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^ I Abstract 

In this work the degree of entanglement of two-photon states gener- 
ated by parametric-down conversion with a phase mismatch is computed 
in terms of the temperature of the nonhnear medium where the interaction 
takes place. The minimum interaction time, named birth time of entan- 
glement, needed before entanglement is attained, is computed in terms of 
the temperature and of the phase mismatch. Our results show that, for 
each value of the phase mismatch, there is a critical temperature above 
which entanglement disappears. We are able to identify the regions of 
parameters, temperature and mismatch, where the states are entangled 
and where they are separable. We also propose a entanglement witness 
based in a dynamical invariant of the nonlinear interaction. 

PACS numbcr(s) 03.67.a, 42.50-p., 42.65.Lm 
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1 Introduction 



A significant number of Quantum Computation and Quantum Information pro- 
tocols use the quadrature amplitudes of the quantized electromagnetic field [1] , 
[5J [3]. The essence of all these protocols rely in the possibility of producing 
entangled photons in a nonlinear medium 0]. 

The Parametric Down Conversion (PDC) [SJIB] is one of the most often used 
processes for production of entangled pairs of photons. However, as for many 
nonlinear optical processes, PDC is ruled by phase matching whose importance 
was early recognized [7] , Theoretical treatment of a phase mismatch that occurs 
during the propagation in the nonlinear medium is complicated and several ap- 
proximations have been made [Sj El [10] , in order to relate the phase mismatch, 
with the amount of quantum correlations in the spontaneous emission of pho- 
tons. 

At room temperatures the average number of photons at optical frequen- 
cies is very small (on the order of 10"'*'^). However, in the microwave part of 
the spectrum the number of thermal photons initially present in the nonlinear 
medium can not be ignored and we expect that they will play a determinant role 
in the amount of entanglement that can be extracted from the interaction. To 
our knowledge a comprehensive study of PDC processes with a phase mismatch 
and at finite temperature of the nonlinear medium has not been addressed. 

In this paper we describe the PDC process with a phase mismatch with a 
simple model that is the quantum transcription of the classical PDC approach 
derived by Blocmbergen |11| . This model is valid when the interaction is weak, 
and when the complex amplitudes of the interacting waves vary slowly with the 
propagating distance in the nonlinear medium. Our purpose is to relate the joint 
effect of the temperature of the nonlinear medium and of the phase mismatch, 
with the degree of entanglement of the down-converted photon states. 

To quantify the degree of entanglement in terms of the phase mismatch and 
of the temperature we use the logarithmic negativity [12]. We show that for 
finite temperatures there is a minimum time of interaction needed, the Birth 
Time of Entanglement (BTE), before the system starts to be entangled and 
that it varies with the temperature and with the phase mismatch. For a given 
phase mismatch, we compute the degree of entanglement with the temperature 
in order to obtain the critical temperature, above which the quantum state 
becomes separable. This study allows us to identify what are the regions of 
joint values of temperature and phase mismatch, where the states of the system 
are entangled and where they are separable. Finally, we propose a feasible 
entanglement witness that is directly connected with a dynamical invariant of 
the system. 

The plan of the paper is the following: In Sec. II we introduce the Hamil- 
tonean of the system with non zero phase matching, and obtain the time depen- 
dent bosonic operators. In Sec. HI we compute the birth time of entanglement 
and the degree of entanglement in terms of the phase mismatch and of the 
temperature and, an entanglement witness is proposed. Finally in Sec. IV we 
present the conclusions. 



2 



2 System dynamics 



The parametric down conversion is a process where an intense pump wave of fre- 
quency Wo, propagates in a nonlinear medium and gives rise to two other modes 
of smaller frequencies lui (the signal) and UJ2 (the idler) obeying the energy 
conservation rule loq = loi + i02. When the interaction is weak, and when the 
complex amplitudes of the interacting waves vary slowly with the propagating 
distance so that they can be assumed approximately constant within a distance 
of a wavelength, it is possible to use the slowly varying envelope approximation 
derived by Bloembergen in a classical context, to describe the nonlinear 
interaction with a phase mismatch S. The quantum version associated to the 
classical Hamiltonean describing this process, is 

H = noJia\ai + huJ2ala2 - %(e-*[("°+*')*+'^lal4 + h.c.) (1) 

where the intense pump mode is treated classically as a coherent undepleted field 
of complex amplitude ao = |ao|e^*[('^°+'')*+'^l and an arbitrary phase if. Modes 
1 and 2 are described by the bosonic operators di and 0,2- The interaction time 
t is the propagation time in the nonlinear medium and the coupling constant 
g is proportional to the second order susceptibility of the medium and to the 
modulus |ao| of the pump wave. 

Solving the Heisenberg equations of motion of the system (see the Appendix) , 
we obtain the time dependent bosonic operators 

ai(T) = e-*™i(^) \a{xt) dio + B{xt) (2) 



B*{xt) aw+A*{xT) aU (3) 



where t = gt is an effective dimensionless interaction time, y ~ S/g is the 
dimensionless mismatch parameter with values in the interval < J/ < 1, 
^'j{y) = '^j / S + y is the dimensionless effective frequency of oscillation of mode 
j and a; = (1 — y^)^/^. The subscript zero refers to initial time. The coefficients 
A{xt) and B{xt) are defined by 

A(xt) = cosh(a;T) H sinh(a::T) (4) 

X 

B{xt) ~ — sinh(a;T) (5) 

X 

Equations ^ and ([3]) extend the well known result derived in [5] to take 
into account the influence of the phase mismatch, and reduces to it when y = 0. 
Because the pump amplitude |ao|, is treated as constant, the solution to the 
parametric interaction ceases to be valid when the number of generated pair of 
photons is such that determines an appreciable depletion of the pump mode. 

The argument of the hyperbolic functions is the squeezing parameter, 

r{xT) = Tx (6) 
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It depends not only on the dimensionless interaction time r, but also on the 
mismatch. Making 5 = in the last equation we obtain the usual squeezing 
parameter r = gt For a given interaction time the squeezing parameter 
decreases with increasing mismatch and we expect that the efficiency of the 
correlated pair production also decreases with y. 

3 Degree of entanglement 

In this work, we are interested in the amount of entanglement we can extract 
from the correlated pair production in terms of the phase mismatch and at finite 
temperature T. 

In the spontaneous PDC emission the signal 1 and the idler 2 modes are 
initially in the vacuum state |0)jo (j = 1,2). The initial state of the composite 
system is the pure and separable state given by the tensor product 1-0)0 = 

|0)lo|0)20. 

When the nonlinear medium is in thermal equilibrium at temperature T, the 
initial state of the two modes is a separable mixed state given by the density 
operator po = pw (8i P20, where pjo = X]^^=o e"^^"^ l'^i)('^il is the thermal field 
density operator of mode j, with /3j = Tiujj/{kBT) and ks is the Boltzmann 
constant. 

The initial mean number of photons in mode j is ftj^iT) = l/(e'^'~^). At 
room temperature and in the optical part of the electromagnetic spectrum 
TijQ << 1. however, in the microwave part of the spectrum fij^ » 1 and 
we cannot ignore the presence of thermal photons in the nonlinear medium. 
The initial average number of photons in the limit T — > is n^o = 0, and we 
recover the vacuum state. 

Both vacuum and thermal states are Gaussian states and since the Hamilto- 
nian H is bilinear in the field modes, the overall output state is also Gaussian and 
the separability criterion [131 HI] is completely characterized by the first and the 

second statistical moments of the quadrature phase operators qj ~ \J'^{0'j'^'^\) 

and pj = --j=={aj — d^), (j = 1, 2). These operators can be grouped together 

in a vector of operators X = {qi,pi, q^iPi)^ ■ 

The first and the second statistical moments of the quadrature phase opera- 
tors, will be denoted, respectively, by the vector X = ((aii), (2:2), {pi), {P2)) and 
by the covariance matrix (CM) a, of elements aij = ■^{xiXj + XjXi) — {xi){xj). 
The first moments can be adjusted to zero by local unitary operations which 
leave invariant entropy and entanglement. The time dependent covariance ma- 
trix cr(T) can be expressed in terms of the three (2 x 2) block matrices air), 
f3{t) and -f{t) 

"<^' - ( iv^ ) 

The diagonal blocks ^(t), /3(r) are the local CM of modes 1 and 2, respectively. 
The ofF-diagonal block, 7(t), encode the intermodal correlations (quantum and 
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classical) between the two modes. 

Hamiltonians that are quadratic in the bosonic operators, possess universal 
quantum invariants, i.e., certain combinations of variances which are conserved 
in time independently of the concrete form of coefficients of the Hamiltonian 
[T5] . These invariants exist due to the symplectic structure of the transforma- 
tion relating initial and time-dependent values of the quadrature components 
operators. Concerning the statistical moments, these universal invariants are 

m 

Ii = det cr(T) ; I2 = det a(r) + det P{t) + 2 det 7(t) (8) 
Defining the quantities 

5(r) =5o + i(det7(T)-|det7(T)|) ; 50=11-^12 + ^ (9) 

the entanglement criterion can be written in a simple form 

S{t) < (10) 

For both vacuum and thermal initial states, we prove in the Appendix, that 
S{t) = So — |det7(t)|, with So ~ "■io"-2o(^io + 1)(^20 + !)• For the initial 
vacuum state, nio = 71,20 ~ and So ~ 0. 

3.1 Birth time of entanglement 

Using the entanglement criterion (|10p . we conclude that the two modes are 
entangled for interaction times tj obeying the inequality 

|dct7(rO|>5o (11) 

where |dct7(ri)| is an increasing function of r. The instant T£, such that 
I det 7 (rg) I = So is the minimum interaction time needed before entanglement 
appears, for obvious reasons it will be named the birth time of entanglement 
(BTE). 

For the vacuum initial state, r > S{t) > 0, therefore the two modes 
are entangled since the very beginning of the interaction, this is, rg = for any 
value of phase mismatch y. On the contrary, for an initial thermal state, the 
condition ([TT|) is attained only after a finite time of interaction, rg, has passed. 

The birth time of entanglement depends on the temperature T, through the 
initial average number of photons nio{T) and h2o(T). To compute the initial 
average number of photons we assume two typical frequencies in the microwave 
region of the eletromagnctic spectrum: i^i = 3.12 x lO^^Hz and 1/2 = ^vi and 
for the coupling constant we assume g = tt x lO^^vi. 

The dependence of the dimensionless BTE, Tg, over the temperature is dis- 
played in Fig.l, for the dimensionless frequencies wi = 2Tn'i/g = 200, W2 = 
2'KV2/g = 400 and for different values of the dimensionless mismatch parameter 
y: y ^ (triangles), y = 0.5 (circles) and y = 0.9 (stars). We observe that, for a 
given mismatch parameter y, the BTE increases with the temperature although 
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this growth slows down at higher temperatures. As we will show, only when the 
average number of correlated photon pairs is higher than a threshold defined 
by the uncorrelated initial thermal photons is the break- even for entanglement 
attained. 

5 - 
4 - 
3 - 

* 

2 - • 
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Figure 1: Dimensionless BTE, T£, over the temperature, for the dimensionless 
frequencies wi = ^il g ~ 200, u)2 = ^2/9 = 400 and for different values of the 
dimensionless mismatch parameter y: y — (triangles), y — 0.5 (circles) and 
y = 0.9 (stars). 

In Fig. 2, the BTE is displayed as a function of the dimensionless mismatch 
parameter y, for different temperatures: T — 50°K (triangles), T = 150°K 
(crosses) and T = 300° K (circles). We observe that, for T > 0, the BTE 
increases with y, and the growth rate is higher at higher temperatures. This 
is not surprising since the squeezing parameter, responsible for the efficiency of 
the photon pair production, decreases with increasing mismatch introducing a 
randomness in the interaction. 

3.2 Logarithmic negativity 

The logarithmic negativity £j\f [121 117] provides a proper quantification of en- 
tanglement in two-mode Gaussian states and is a computable measure of entan- 
glement. 

In Fig. 3 the degree of entanglement £j\f is displayed in terms of y for Ti = 
4.543 and for T = 0°K (crosses), T = 50°/\ (triangles) and T = 300°K (circles). 
The same dimensionless frequencies as in Fig. 2 are assumed. We observe that, 
for a given temperature, the degree of entanglement decreases with increasing 
y. Moreover, for T = 300°K and for y > 0.9, = 0, this is, the two modes 
are in a separable state though for T = 50° if the system is still entangled for 
y = 0.99. 
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Figure 2: BTE is displayed as a function of the dimensionless mismatch param- 
eter y, for different temperatures: T = 50° K (triangles), T ~ 150° K (crosses) 
and T = 300°K (circles). 



As a matter of fact, when the initial temperature of the nonhnear medium 
rises, there are initially more uncorrelated thermal photons in the modes. We 
expect intuitively that, for a given interaction time, any entanglement should 
vanish at some critical temperature Tc- We have computed £j\j-, for ~ 2.881, 
in terms of the temperature and for different values of the mismatch parameter 
y, in order to find the critical temperature for each y. The results are displayed 
in Fig. 4 where the points {y, Tc) are marked with circles. The solid curve was 
obtained by interpolation. The quantum states of the two modes characterized 
by the parameters y and T in the shaded region under the curve are entangled 
states. On the contrary, above the curve these parameters define separable 
states. We observe that there is a kind of trade-of between the phase mismatch 
and the critical temperature: greater is the phase mismatch smaller is the critical 
temperature and vice versa. 

Another interesting conclusion is that it is possible to have entangled states 
at room temperature for small enough phase mismatch. For instance, in the 
case considered in the Fig. 4, the modes are entangled at T = 300°K for S < 0.3. 

3.3 Entanglement witness 

In Fig. 5 is displayed the average number of photons {n{T)) ~ ((ni(T)) + 
(n2(r)))/2 produced in the PCD process after a time of interaction = 2.881, 
in terms of the temperature T and for different values of the mismatch param- 
eter y = (dashed-dot line), y — 0.5 (dashed line) and y = 0.7 (dotted line). It 
is also displayed the invariant \/Sq (solid line) in terms of the temperature. The 
first three lines intersect \/So at different values of the temperature T. Com- 
puting the temperature at the intersection points we find out that they coincide 
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Figure 3: Degree of entanglement Ej^ , in terms of y for = 4.543 and for 
T = Q° K (crosses), T = bQ°K (triangles) and T = 300° A' (circles). The same 
dimensionless frequencies as in Fig. 2 are assumed. 



with the corresponding critical temperature Tc for each y. For temperatures 
such that the average number of photons is in the shaded region of the figure, 
the system is entangled. Under the solid curve the average number of correlated 
photons produced in the interaction arc not enough to attain the break- even, 
given by the square root of the invariant 5*0, this is, they arc not enough to 
exceed the uncorrclatcd photons present in the initial thermal state, and the 
system is separable. 

We can then define a simple entanglement witness 

W{T) = v/Sio(T)n2o(r)(nio(T) + l)(n2o(r) + 1) - (n(T)) (12) 

As usual, the negative value of yV(r) indicates entanglement. To check wether 
the two modes are entangled, it will be enough to measure the average number 
(n(T)) of correlated photons generated in the two modes and the equilibrium 
temperature T of the nonlinear medium, since this define the initial average 
number of photons UjolT), (j = 1,2) in each mode. 

4 Summary and concluding remarks 

We have computed the joint effect of a phase mismatch and of the temperature 
in the degree of entanglement in PDC processes described by Hamiltonean ([T|). 

Using the time invariants associated to that hamiltonean we computed the 
minimum interaction time needed before entanglement appears and have shown 
how it behaves with the phase mismatch and temperature. 

We have also shown that for each phase mismatch there is a critical tem- 
perature, above which the entanglement disapears. This study enables us to 
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Separable 




Figure 4: Critical temperature Tc for different values of y (circles). The solid 
curve was obtained by interpolation. The shaded region under the curve corre- 
spond to values of parameters (T, y) where the states are entangled. Above the 
curve these parameters define separable states. 



draw the border line between separable and entangled states in the plane of the 
parameters, phase mismatch and temperature. Depending on the interaction 
time, and on the phase mismatch it is possible to obtain entanglement at room 
temperature. Finally we propose a feasible entanglement witness. 



5 Appendix 

The Heisenberg equations of motion for modes 1 and 2 are 



dt V a 



fli 



it 



ige 



-i[2{u+5)t+ip\ 
iUJ2 



(13) 



jointly with their hcrmitian conjugates, which are decoupled from these two. In 
order to obtain an autonomous system of equations we define new annihilation 
operators bi and &2, by the canonical transformation 



bi 
bl 



gi(w+i5)t g 

Q g-i(a;+(5)t 



al 



(14) 



Physically this corresponds to move to a reference frame that rotates with the 
frequency {ur + 6). The Heisenberg equations for the new operators bi and b^ 
are then given by 



dt 



i(u! — LOl + 5) 



ige 



-lip 



-ige 



lip 



i{uj ^ uji — S) 



(15) 
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Figure 5: Average number of photons {n{T)) = {{ni{T)) + (7i2(T')))/2 produced 
in the PCD process after a time of interaction ti = 2.881, in terms of the 
temperature T for different values of the mismatch parameter: y — ( dashed- 
dot line), y ~ 0.5 (dashed line) and y ~ 0.7 (dotted line). The invariant ^/Sq 
(solid line) is displayed in terms of the temperature T. 



They are linear time independent differential equations and can be readily in- 
tegrated. The only solution of eq. lfTS)) with physical interest is given by eqs.([2l 
[3]) of Section 2, and corresponds to the condition S < g where the bandwidth is 
small compared with the coupling parameter. 
It is immediate to prove that 



{hi{t)) - {n2{t)) = nw - n2Q 



(16) 



is an integral of motion. This is a generalization, for nonzero mismatch, of the 
analog integral of motion derived in [5]. 

Using the time dependent bosonic operators given in eqs.© and ([3]) we 
compute the entries of the covariance matrix a for the thermal initial state 



M = (n,(r)) + i ; 0- = l,2) 



(17) 



(Ti2(t) = (721 (t) = crsiir) = (T43(t) = ; (j = 1, 2) 
0-24 (r) = 0-13(7-) ; CT23(t) cri4(T) 



CTi3(t) = (»^10 + »^20 + 1) 



cti4(t) = (nio + ^20 + 1) 



S{xt) 

X 

S{xt) 

X 



sm{eT)C{xT) 



cos(6'r)y 



S{xt) 



cos{0t)C{xt) + 5^^^i^5(xT) 



where 6 = {uj[ + uj^) and w' = ujj + y, j = 1, 2. 



(18) 
(19) 

(20) 
(21) 
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The CM for the vacuum state is easily obtained, by making nio = 7120 = 
in the above equations. 
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